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Massless two-loop self-energy diagram: 
Historical reviewo 

A. G. Grozin 

Abstract 

This class of diagrams has numerous applications. Many interesting results have 
been obtained for it. 



1 Introduction 

We consider the integral (Fig. [1]) 
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(1.1) 

in d = 4 — 2e-dimensional Euclidean momentum space. It has a long and interesting 
history. For many years, most of the information we had about perturbative quantum 
field theory was coming (directly or indirectly) from this integral. All massless three- 
loop self-energy integrals (with integer indices) reduce to 6 master integrals PQ, 5 
of which are particular cases of / (jl.ip . Only one master integral (the non-planar 
one) does not reduce to /; however, using the gluing method pQ, one can easily show 
that its value at e = is equal to the ladder integral, which reduces to 1(1, 1, e, 1, 1). 
At four loops [2j, 15 master integrals (of 28) reduce to /, and thus can be easily 
expanded in e up to high powers. 





Figure 1: Two-loop self-energy diagram in momentum and coordinate space. 



* Extended version of the talk at the international conference Advances of quantum field theory, Dubna, 
October 4-7, 2011. 
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This integral in coordinate space (Fig. [I]) 

d d Xl d d x 2 



/(ai,a 2 ,a 2 ,a 4 



[(Xl - X) 2 ]" 2 [(X1 - X 2 ) 2 P[(X 2 - X) 2 f^{xl)^ 

I(a 2 , 0,4,0,3,0,5,(11) , (1.2) 



has the same form (jl.ip if we rename Xj —tpf, here trivial T-functions from Fourier 
transforms are not explicitly shown, and 

ai = ~-<k- (1-3) 
We can perform inversion of the integration momenta in (jl.ip 

fr. - _*± J- 2 - J_ _ d k 'i (] ? . ] r .\2 _ ( k l Z ^ 

K * ~ JJ2 ' K i - U2 ' " Kl - /J./2\d ' K2 > ~ U2 y2 ' 

and obtain 




Inversion relations can be also derived in coordinate space, of course. 
The integrals 




where dashed lines have integer indices, have been calculated in [3] via Gegenbauer 
polynomials. Of course, now we know that it is trivial to calculate them using IBP 
(Seet.E). 



2 Integration by parts 

The IBP relations for this particular class of integrals first appeared in [4j. They 
are described in the text below the formula (15); this formula is the homogeneity 
relation (which is a consequence of the IBP relations). Soon IBP relations evolved 
into a fantastically universal and efficient method for reducing all scalar integrals of 
a given topology to a few master integrals pQ. 

The IBP relations allow one to trivially reduce integrals with integer indices in 
the left (or right) triangle to one-loop integrals expressible via T-functions: 
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If 03 is not integer, things are more difficult. The combination [lj of the IBP 
relations 



[(d - 2a 3 - 4)3+ + 2(d - a 3 - 3) 



a 3 )— = 21+(5~-2-3 + )— < 



03 , 



(2.2) 



allows one to shift 03 by ±1 (if all integer indices are 1, all integrals in the right-hand 
side of (|2.2p are trivial). 



3 Uniqueness 

Many interesting results for massless self-energy integrals were obtained using the 
method of uniqueness [U G2 O [7] (see also the textbook [$]). It is based on the 
following relations. In coordinate space 




a% + ci2 



ai 



02 



a\ + a% 



(3.1) 



02 



(in momentum space the second formula becomes trivial, and the first one contains 
some factor; these combinations of T- functions from Fourier transforms are not ex- 
plicitly shown here). The main element of the method is the star-triangle relation 
which is valid if a\ + 0,2 + (13 = d: 



(3.2) 



where 04 are defined by (jl.3p (note that a\ + 02 + a 3 = d/2). It can be easily derived 
using inversion. 

Kazakov [7] has calculated a non-trivial integral 




1(a) 



(3.3) 



(the dashed lines have indices 1) via hypergeometric functions of the argument —10- 
It has a symmetry property 



1(1 + a) = J(l-o-3e) 
(see Sect.HD; 1(0) is known. The IBP relation ([22]) gives 



1(1 + a) 



l-a-2e 
a + e 



1(a) 



(l-2a- 3e)r 2 (l - e)T(-a - e)T(a + 2e) 
(a + e)r(l + a)r(2-a-3e) 



(3.4) 



(3.5) 



1 Earlier a particular case of this family 1(e) has been calculated via hypergeometric functions of 1 [S], 
see Appendix O 
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If we write 1(1 + a) via a new function G(l + a) 

/(l + a)=2 r2 ' 1 - e > r '-°- e ' r '° + 2e ' 
v ; r(i + o)r(l - o - 3e) 

then this recurrence relation becomes simpler: 



G(l + a) 



G{1 + a) 



-Ma) + 



1 



1 



+ 



1 



1 — a — 3e a — l + 3e \ a + e a — 1 + 2e 

Writing this function as a sum over its poles 



G(l + a) = J2A 1) 



1 



+ 



1 



n + a + e n — a — 2e 



+£/■?> 



1 



+ 



1 



n + a n — a — 3e 



n=l \ / n=l 

(where the symmetry (|3.4p is taken into account), we obtain recurrence relations for 
the residues: 

n + e 



f(i) 

J n 



n + l-2e 



f W f(2) 
J n+l ' Jn 
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n + 1 — 3e 



f (2) 



Their solution is 



/ » 1> " ( - 1) " n ^T7r <!l(E) ■ / » 2) - ( - 1) " r( "r ( l)" 3£> <:2(£) ' 

where the constants are obtained from the initial condition: 

T(e) r(e)r(l - e)r(l + e) 



ci(e) 



c 2 (e 



T(2 - 2e) 
Therefore we arrive at 

r 2 (l - e)r(e)r(a + 2e)r(-a - e) 



r(2 - 2e)r(l - 2e)r(l + 2e) ' 



1(1 + a) = 2- 



EC- 1 )' 



r(2 - 2e)r(l + a)r(l - a - 3e) 
>r(n + l-2e) 



.n=l 



r(n + e) 
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+ 



n + a + e n — a — 2e 



r(l-e)r(l+e) 
r(l - 2e)r(l + 2e) 



D-d 



n r(ra+l-3e) / 1 



n=l 



r(n) 



+ 



n + a n — a — 3e 



This result can be written via hyper geometric functions: 



1(1+ a) = 2 
1 



r(e)r 2 (l - e)T(a + 2e)T(-a - e) 
r(2 - 2e) 

l,2-2e,l + o + e 



T(2 - 2s) 



T(l + e)r(l + a)T(l - a - 3e) 



1 + a + e 
1 



3-^2 



+ 



1 - a - 2e 



3^2 



L + e,2 + a + e ' 

l,2-2e, 1 -a-2e 
l + e,2-a-2e 





— cos(7re) , 







(3.6) 



Kazakov [BJ [7] also derived several terms of expansion of I(a,i) with ai = 1 + nje 
in e using symmetry properties of /; we shall discuss this expansion in Sect. HI 
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4 Symmetry 



Symmetries of the integrals (jl.ip which follow from inversion (|1.4p . duality between 
x and p space (|1.2p , and the star-triangle relation (|3.2p were considered in [3] 

Gorishnii and Isaev |10j discovered the tetrahedron symmetry group S4 of the 
integrals /. Let's consider the vacuum diagram in Fig. [21 all lines have mass m. If 
we integrate in the momentum of the line 6 last, then 

1 f F(k 2 )d d k 



n d/2 J (k 2 +m 2 ] 



where F(k 2 ) is the self-energy diagram with external momentum k obtained by cut- 
ting the line 6. Its asymptotics is 

^ fl 2 \ I(ai,a 2 ,a 3 ,a^,a 5 ) 
F(k — > 00) 



^2 s ja 1 +a 2 +a 3 +a4+a 5 -d 



(it comes from the hard region; other regions give contributions with different powers 
of k 2 ). Hence the vacuum diagram has the ultraviolet pole 

1 1(01,02,03,04,05) 



r(d/2) 01 + a 2 + a 3 + a 4 + a 5 + a 6 - §d 



(other regions produce poles at different places). But we can equally well cut some 
other line, I\jy must remain intact. Therefore I{a\, 02, 03, 04, 05) (which also depends 
on d) can be considered as a function of 01,02,03,04,05,06, where a& is defined by 

3 

«i + a 2 + a 3 + o 4 + a 5 + a 6 = -d , (4.1) 
with the tetrahedron symmetry. 




Figure 2: The tetrahedron diagram. 

Gorishnii and Isaev also considered symmetry relations following from the star- 
triangle relation (13. 2p (which were discussed in [HE]). Taken together, these sym- 
metry transformations are sufficient for generating the complete symmetry group of 
the integrals I. But they could not identify this group. A complete solution of this 
problem was obtained in [12"] . 
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Let's introduce notation 



I(a 1 ,a 2 ,a 3 ,a 4 ,a 5 ) 



10 



net 

i=l 



1/2 



/(ai,a 2 ,a3,Q4,Q5 ; Q6) 

(d - 3)r2 (| - 1) 



(4.2) 



where a§ is defined by (|4.ip , aj by (jl.3p , indices of the vertices (Fig. [2]) are 

d d d d 

a-i = ai+a 5 +a 6 -- , a 8 = a2+a 4 +a 6 -- > 09 = 03+04+05--, a 10 = a 1 +a 2 +a 3 -- , 



and 



G(a) 



r(q) 
r(o) 



The pre- factor in (j4.2j) is chosen in such a way that T-function factors in (|3.ip . (|3.2p 
cancel in symmetry relations. 

The full symmetry group is generated by 3 transformations. The first two gener- 
ators are elements of the tetrahedron group: 




The last one comes from uniqueness. First we introduce an extra dot on line 3 to 
make vertex 10 unique; then use the star-triangle relation (13.2|) : and then combine 
two lines: 



6 6 




6 





I. e., the third generator is 
10 





(4.5) 



The structure of the group becomes apparent if we introduce new variables: 
I(a 1: a 2 , a 3 , a 4 , a 5 , a 6 ) = I{b u b 2 , b 3 , b 4 , b 5 , b 6 ) , 



(4.6) 
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Then our 3 generators transform i(b\, b 2 ,b 3 , 64, 65, b§) to 

Pil(b 1 ,b 2 ,b 3 ,b i ,b 5 ,b 6 ) = I(bi,be,b2,h,h,h) , 
P 2 I(bi,b 2 ,b 3 ,b4,b 5 ,b 6 ) = I(b 3 ,h,bi,b6,h,h) , 
P 3 I(b 1 ,b 2 ,b 3 ,b i ,b 5 ,b 6 ) = I(b 3 ,b 2 ,bi,b A ,b 5 ,b & ) . 

We can combine them into 3 better generators 

Qi = (P 3 P!) 2 P 3 P 2 Pi , Q 2 = PfP 3 Pi , Q 3 = (P3P2P?) 2 ; 

they transform I(bi,b 2 , b 3 , 64, 65, be) to 

Qil(bi,b 2 ,b 3l b4,b 5 ,b 6 ) = I(b 2 , b 3 , 64, 65, 6 6 , 61) , 
Q 2 I(bi,b 2 ,b 3 ,b4,b 5 ,b 6 ) = I(b 2 , b%, b 3 , 64, b 5 , & 6 ) , 
Q3l(h,b 2 ,b 3 ,b i ,b 5 ,b 6 ) = I(bi,b 2 ,b 3 ,h,b5,h) ■ 

The first two generate the symmetric group Sq of permutations of 6 variables &i,...,6; 
and Q 3 generates Z 2 . So, the symmetry group of the integrals / is Sq x Z 2 [12]; it 
contains 6! • 2 = 1440 elements [TT] , 
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The most useful information is the expansion of 1(01,02,03,04,05) (jl.ip in e at 
d = 4 — 2e and ai,...^ = 1 + 0(e). This means that all Oj are 1 + 0(e); the same is 
true for &, (note that X^=i &i = 1^)- The function (|4.6p is invariant with respect to 
S"6 x Z2; therefore, the expansion can be written entirely via invariants of this group. 
The invariants are 



1 



d 



i=i 



and 



i=i 



for n = 2, 3, 4, 5, 6. The total degree of I3 and 1$ must be even, because they change 
their signs under Q5: 



1(01,02,03,04,05,06) 



jtl p2 r«3 r«4 r«5 pa 
«2«3M«5*6-'l J 2 J 3 J 4 J 5 J 6 • 



(4.7) 



i3+*5 even 



All unknown coefficients C^^isis needed for expanding up to e 4 can be fixed by 
considering the integrals (|2.ip (the indices in the left triangle are 1). IBP reduces 
these integrals to T-functions, and hence all these coefficients are expressed via £ n . 
So, the symmetry allows one to obtain, practically for free [6| 17] [TT| H2]. 



5 

1(01,02,03, a 4 , 05, a 6 ) = 6C3 + I8C4/1 + 3(5 ( if + 7^2 



15(7C 6 + 2^)1? + 3 ( — Ce - C3 ) hh 



25 



439 



211 



9 ( — Cr + 2OC4C3 Mi + 3 ( — Ct - 6C4C3 ) I{ h + -(7 [ -jli ~ 7/4 ) + 



35 



(4.8) 



5 Gegenbauer polynomials 



If indices of 3 lines forming a triangle are 1, I reduces to T-functions by IBP (|2.ip . 
In a more general case when 2 adjacent lines have indices 1 I can be expressed via 
hypergeometric functions of 1. Due to the tetrahedron symmetry, it does not matter 
which 2 adjacent lines have indices 1 (Fig. [3]). These integrals were calculated by 
Kotikov [13] by an ingenious use of x-space Gegenbauer polynomials [3]. 






Figure 3: Integrals with 2 adjacent lines having indices 1. 
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Let's consider the first of these diagrams in x-space: 



I(c,b, 1,1, a) 



G 2 (l)G(a)G(b)G(c) 
G{a + b + c + 2- d) 



A(a,b,c) , 



(5.1) 



A(a, b, c) 



1 



d d x d d y 



v d J {y 2 ) a [{z-y) 2 ] x [{z-x) 2 ] b {x 2 ) c [{x-y) 2 } x ' 
where A = d/2 — 1 and z 2 = 1. Expanding the propagator as 



1 



1 ^T(A + n) 



[(y r(A)^ o m 



E 



-v 



fll. ..fir 



fll...fl„ 



(where yW— j s the traceless part of y^ 1 ■ ■ ■ y^ 1 ™) we can integrate in d d y: 
1 ^ 2 ™r(n + A)^ 1 ... M , 



(5.2) 



A(a, b, c) 



■E 



T 2 (A)(a-l)^ n\ ir d / 2 J {x 2 Y[{z-xy} b 



X X' 



Ill-.-tln 



n + A — a + 1 



9(1 -x 2 ) + fl(x 2 - 1) 



1 



(x 2 ) n+A 7 n + A + a-1 
Now the integral in d d x can be calculated 

1 



6{l-x 2 ) + 



9(x 2 - 1) 



A(a, 5, c) 



r(A)r(2A)r(b)r(6- A)(a-1) 



^ 2 n r(n + A) ^ r(m + n + b)T(m + b - A) 



n=0 



1 



1 



+ 



1 



m=0 



mlT(m + n + A + 1) 



n + \ — a + 1 Vm + n — a — c + A + 2 m + n + 6 + c— 1 



1 



1 



+ 



1 



n + A + a — 1 \m + n — c + A + 1 m + n + a + 6 + c — 2 

It appears to be possible to transform this result into a form containing only 
single sums [15] : 

A(a,fe,c)= ,/\T^ 2 ^ (5.3) 



r(f-l) (a-1) 



where 



An 



&)r i 



r(f- a 



r(6) [d-2a 

c+l)r(a + 6 + c-f-l) 



(5.4) 
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| — a + 1, d — a — b — c+1 
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r(c+f-i)r(f-6- c + i) 
r(|- c )r(6 + c -|) 



3-^2 



1 



2a + d - 4 



T(c)T(d-b-c) 



3F2 



d-2, f-a,6 + c-l 
|-a + l,c + |-l 

d-2,a+|-2,f-e 
a + I — l,d — 6 — c 
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+ 



A, 



r(2 



c)Y{a + b + c-2) 



r(a + c +|-2)r(| 



:+2) 



3-^2 



2,a+f 



2,a + 6 + c-2 



a + 



l,a + c+ | - 2 



r(i-b)r(i-c)r(^-a)r(a + |-2)r(^-b)r(a + b + c-|-i) 
r(d - 2)r(o + c - i)r (a + 6 - 1) r (f - a - b + 1) r (| - 6 

2r(i - 6) 



+1) 



(2a + d - 4)r (6 - f + 1) 
T(2-a-c)r(a + 6 + c- 



1 



r(3-a-6-c)r(a + c+- 

r(i-c)r(6 + c -|) 



+ 



r(c-f + i)r(d-6-c 



-2) 

3-^2 



3-^2 



-2,a + f 
a + ^- 



2,a + 6 + c-2 



i-2,a + |- 
a+ | - l,d 



l,o + c+f 
1 



o a, 

Z, 2 

-6-c 



(5.5) 
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r(l-6)r(2-a- c)r(f-a)r(a + f-2)r(f-6)r(fe + c-f) 

r( c )r(d - 2)r (6 - 
2r(i - b) 

(d - 2a)T (b - | + 1) 



-f + 2)r(a-6+i-l)r(|-a 

r(i- c )r (6 + c-j) / d- 



r(2-6-c)r( c + f -l) 



r(2 



c)r(a + 6 + c-|-l) 



3-^2 
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| -a + l,d 



a, 5 



c+2) 

- a, 6 + c 
d 
2 

C+l 



9 ^ 
Z ' 2 

|-a+l,c+|-l 



r (a + c - f) T(d - a - b - c + 1) 

r(i - 6)r(i - c )r(2 - a - c)r (g - a) r (a + 1 - 2) r ( 



b-c + l 



(5.6) 



6) 



r(i-o)r(o)r(d-2)r( 



+ 



2r(l - b) 



T(b- 
3F2 



+ 1) 

d- 



r(2 



a-6 + l)r(f 
c)r (a + & + c - 



+ 2) 



1 



(d 

9 <i 
Z ' 2 



a - 

u, 2 



I - a + l,d 



r(i-c)r(6 + c -|) 



2a)r(a + c- f)r(d 

a — c + 1 
a — 6 — c + 1 

d 1 

3^2 



6 -c+l) 



Z-2,a + f - 
a + f - l,d 



J > 2 
6- 



(2a + d - 4)r (c - I + 1) r(d - 6 - c) 
r(i - 6)r (| - a) r (a + 1 - 2) r (f - 6) r (6 + c - 1) r (a + 6 + c 



(5.7) 



1-1) 



+ 



2r(i - 6) 
r(6-| + i) 

r(2 - a 



r(i - o)r(o)r(c)r(d - 2)r ( c + § - 2) 
r(i -c 



,)T{b + c-l) 



(d - 2a)r(2 - b - c)r (c + I - 1) 
-c)r(a + 6 + c-|-l) 



3-F2 



d-2,1 
f-a4 



(2a + d - 4)r(3 



3-^2 



d- 2, a + 
a + ^- 



6-c)r(a + c+§ -2) 
-2,a + 6 + c-2 



1, a + c + 



In particular, for the integral 1(a) (|3.3p we obtain 



/(a) = 2r 



i)r( 



l)T(a-d + 3) 



2r(| 



.(d-2o-4)r(a + l)r(§d 



3-^2 



l,d-2,a-f + 2 
a + l,o- I + 3 



- a, 6 + c — 1 
l,c+| - 1 



(5.8) 

(5.9) 

7r cot 7r(d — a 
T(d-2) 
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this form is equivalent to (|3.6p . though mathematical proof is unknown. 
No expression for I with unit indices of 2 non-adjacent lines is known. 



6 Solving IBP for 3 non-integer indices 

Expressions for I with 2 indices of adjacent lines equal to 1 via hypergeometric 
functions of 1 were also derived by guessing the solution of IBP for such integrals |14| . 
Let's consider 

v«2 



/(ai,a2,a 3 ,a 4 ) 




«3 



(6.1) 



where a 4 = ai + 02 + o 3 — i ■ The IBP relations are 

a\I{a\ + 1, 02, 03, 04 + 1) — (oi + 02 — d + 2)I(ai, 02, a 3 , 04) 

'aiG(ai + l,a 3 + 1) „, n 
■ - G(a 2 ,a 3 + 1) 



3 G(1,04 + 1), 

V 04 — a 2 + 1 
(a 3 - d+ 2)/(ai,a 2 ,a 3 ,a 4 ) 

(ai + a 3 - d + l)(a 2 + a 3 - d + 1) 



a 3 - d/2 + 1 
o 3 (a 3 + 04 - d + 1)G(1, a 4 ) 



7(ai,a 2 ,a 3 - l,a 4 - 1) 



G(ai,a 3 + 1) G(a 2 ,a 3 + 1 



04 - a 2 



+ 



04 — ai 



(6.2) 



where G(ai,a 2 ) = G(ai)G(a 2 )G(a"i + a 2 ) is the standard massless one-loop self- 
energy. 

If we express J(ai, a 2 , a 3 , 04) as 



d- 3 



-J(ai,a 2 ,a 3 ,a 4 ) 



a 3 a 4 G(l, a 4 + 1) 

= G(ai, a 2 + 1)5(1 - ai — 1, a 2 - 1, f + ai - a 4 — 2, a 4 - a 2 ) + (ai O a 2 ) 
via a new function S^ai, a 2 , a 3 , 04) satisfying 

5(ai, a 2 , a 3 , 04) = 5(a 2 , 01,03, 04) = — 5(a 3 , 04, ai, a 2 ) , 

a/ \ 1 , + a 3)(«l +«4) p , . v 

aii(ai,a 2 ,a 3 ,a 4 ) = H ■ ■ ■ b{ai - l,a 2 ,a 3 ,a 4 ) , 

a\ + a 2 + a 3 + 04 

then (|6.2p holds. The solution of (|6.4p can be written as 

7r cot 7ra 3 1 a 2 + a 3 



(6.3) 



(6.4) 



5(ai,a 2 ,a 3 ,a 4 ) 
where 
H{ai, 02,03,04) 



F(ai+a 3 , -a 2 , -a 3 ,a 2 + a 4 ) (6.5) 
1/(01,02,03,04) a 3 a 2 o 3 

_ r(l + ai)r(l + a 2 )r(l + a 3 )r(l + a 4 )r(l + fll + q 2 + q 3 + 04) 
r(l + ai + a 3 )r(l + at + a 4 )r(l + a 2 + a 3 )r(l + a 2 + 04) 
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and 



F(ai,a 2 ,a 3 ,a A ) = 3 F 2 



Expansion of 



1 + 




1, -ax, -a,2 
1 + a 3 , 1 + a 4 

1 + a 



1. 



(6.6) 



1 + 6) 



to all orders in a, b at e = can be expressed via ( 2 n+i [33] (some particular cases 
were known earlier [JJ [TT| [12] ) . 

Elegant symmetry-based methods to derive several terms of e expansion of hyper- 
geometric functions (|6.6p algebraically are presented in |14J . However, these methods 
cannot be extended to higher orders. They are no longer necessary: the algorithm 
constructed in |15| allow one to expand such hypergeometric functions to any or- 
der, results are expressible via multiple ( values. The knowledge of expansions of 
these 3^2(1) allows one to reconstruct all unknown coefficients in the general expan- 
sion fl47j) up to e 9 [HCE]. The e 5 term 



I(cti, a,2, a 3) 04, a 5 ,a 6 
/378 
l~5~ 



+ 3 I ^Css 



3/9 
2U C53 



33523,. , , \ * /54 A 1009 
4023 , ?rr \ T T 2,J^ 1083 



Cs + 42C5C3J I'lh 
Cs + 8(5(3) hh + 



(6.7) 



is the first term where a depth-2 value £53 appears. 

In particular, we obtain several equivalent results for 1(a): 



(d - 3)(d - 4)I»r (Id- a- 4) 



2I*(g-l)r(a-d + 3)r(g 

3-^2 



1) 



1(a) 



3d -2a -10 _ / l,g-2,a-d + 3 
d — a — 3 
3d - 2a - 10 



d — a — 3 
a - 1 



■ 3-^2 



a, a — d + 4 
1,1 — a, d — a — 3 



2' ' 



1 ) + Air cot 7r(a - d) - 2 
1 . d d - 4 

1 ) + ^7T COt 7Ttj + 

d — a — 3 



(6.8) 
(6.9) 



d - 2a - 2 



3-^2 



l,o- §d + 5,a - f + 1 



r,a 



+ 2 



H+^cotvrf-2-^-^ 



a 



1 



d - 2a - 2 



3-^2 



l,f -2,f -a-1 
|d — a - 4, 1 - a 



1 ) + j4"7T COt 7T 



(6.10) 
(6.11) 



where 



A 



r(o)r (|d - a - 4) 
r(d-4)r(^-i) 
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A curious integral belonging to the current class was considered in |17j . The 
symmetry allows one to write it in 12 equivalent forms: 

1(1, 1, A, 1, A, A) = 7(1, 1, A, A, 1, A) = 7(1, A, 1, 1, A, A) = 7(1, A, 1, A, A, 1) 
= 7(1, A, A, 1, 1, A) = 7(1, A, A, 1, A, 1) = 7(A, 1, 1, A, 1, A) = 7(A, 1, 1, A, A, 1) 
= 7(A, 1, A, 1, 1, A) = 7(A, 1, A, A, 1, 1) = 7(A, A, 1, 1, A, 1) = 7(A, A, 1, A, 1, 1) (6.12) 

where A = d/2 — 1. It reduces to 1(1, A, A, A). However, it cannot be directly cal- 
culated using the above formulas: one of the arguments should be shifted by x, and 
the limit x — > should be taken. In the paper |17j . recurrence relations shifting d by 
±2 were derived (they were used in [IE]). This method became popular later. 

7 Mellin— Barnes representation 

The integral I can be written as (k 2 = 1) 



at? -^4 



where 




I(a l ,a 2 ,a 3 ,a A ,a 5 ) = - m J [( ^ _ fc)2]ai (fc2)a5 V((h - k) ,k,), 

vm - kf, kD = -L ! r „, , , 91 ^ 2W , . (7.i) 



TT d / 2 J [(k 2 - k) 2 ]^[(h - k 2 ) 2 ] a3 (^) ai ' 

Substituting the Mellin-Barnes representation of the one- loop vertex n " 1 



V((k 1 -k) 2 ,k 2 ) = j—^j dz 1 dz 2 [(k 1 -k) 2 }^(k 2 )^v(z 1 ,z 2 ) (7.2) 

(where v(z\,z 2 ) is a combination of T-functions), we can easily calculate the loop 
integral in k\ [2lffi 



1(01,02,03,04,05) = , n _.^ 2 I dzidz 2 G(a\ - zi,a 5 - z 2 )v(zi,z 2 ) (7.3) 



1 

(2ni) 

(G(a\,a 2 ) is the standard massless one-loop self-energy integral). The result is 

1 



1(01,02,03,04,05) 



27rz) 2 r(o2)r(a4)r(a 3 )r(d - a 2 - 04 - a 3 ) 
I - 04 - a 3 - z 1 )Y(i 



. Y{-z x )T{$ - 04 - a 3 - z x )T(i - a x + z x 
x / dz\ dz 2 



T(-z 2 )T(\ _ 02 - 03 - z 2 )Y(\ - a 5 + z 2 ) 



T(a\ - zi) 

'd 
2 



r(a 5 - z 2 ) 

T(a 1 + o 5 - \ - zi - z 2 )T(a 3 + z\ + z 2 )T(a 2 + a 4 + a 3 - f + zi + z 2 ) ^ ^ 
T(d -01-05 + ^1 + 2:2) 



2 For d = 4 and = 1 it was obtained in [20 . 
3 For d = 4 and a, = 1 this was also done in 1201. 
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This double Mellin-Barnes integral can be expressed via double sums [21j . First 
we close the z\ integration contour to the right. There are 3 series of poles: T(—zi), 
r(| — 04 — 03 — z\), T(ai + 05 — I — z\ — Z2), and we obtain 3 sums over residues 

J = h + h + h ■ 

Then we close the Z2 integration contour contour to the right, and get double sums: 

/ = /ll + /l 2 + /l3 

+ I 2 1 + I 2 2 + h 3 

+ ^3 1 + ^3 2 + ^3 3 + 3 4 + ^3 5 • 

These nested sums belong to the classes which can be expanded in e to any 
order in terms of multiple ( values by the algorithms constructed in [15) . These 
algorithms were implemented in the packages NestedSums [22] (in C++ with GiNaC) 
and XSUMMER [23 1 (in FORM). Therefore, expansion of the integrals / to any order in e 
can be written in terms of multiple £ values [21] . 

Acknowledgements. I am grateful to P. A. Baikov, D. J. Broadhurst, K. G. Che- 
tyrkin, A. I. Davydychev, M. Yu. Kalmykov, A. V. Kotikov, V. A. Smirnov for numer- 
ous discussions of various questions related to the present topic; to Yu. M. Pismak, 
A. P. Isaev, R. N. Lee, N. A. Kivel for constructive comments; to T. Huber, D. Maitre 
for their help in using HypExp and HPL; and to D. I. Kazakov and the members of 
the organizing committee for organizing the conference and inviting me to present a 
talk. This work was supported by the BMBF through Grant No. 05H09VKE. 



A Dispersive calculation of 1(e) 



The integral 



, -A- . 
/ \ 
\ / 

' -V- ' 



(Al) 



with all indices equal to 1 was considered in [S]. By dimensionality it is ~ (k 2 



-3e 



— s — iO (s > 0) this factor becomes s ~ 6e e 6me ^ its imaginary part is thus 
s~' ie sin(3-7re). On the other hand, this imaginary part can be calculated via Cutkosky 
rules (Fig. [3J. 



At k 1 

-3e 




Figure 4: Two- and four-particle cuts. 

The two-particle cuts contribution contains the one-loop triangle (with one non- 
integer index) having two legs on-shell; it is expressible via T-functions. The d- 
dimensional two-particle phase space also reduces to T-functions. The integral over 
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the four-particle phase space can be calculated in terms of 3F3 of unit argument, and 
the result is [9] 



1(e) = (d -4) 

2r(|-i) 

(3d - 8)r(2d - 5) 



r(f-i)r(d-3)r(5-fd) 
r(3-j) 

l,d-2, |d-4 
2d-5,|d-3 



(A.2) 



5^2 



1 - T(3 - d)T (2 - |) cos(vrd) 



Using IBP (|2.2p we can easily obtain I(n + e) for all integer n. However, this 
method cannot be generalized to I(n + me) with m/1, unlike (|5.9|) . 



B Expanding hypergeometric functions in e 

As we have seen, a large class of integrals I can be expressed via hypergeometric 
functions of the argument 1 (or —1) with indices which are integer at e — > 0. They can 
be expanded in e using the algorithm A from [15] ; there are also other algorithms [2i\ 
25j. It is easy to expand them up to quantities of transcendentality level 8 using the 
Mathematica package HypExp [24J (which uses the package HPL |26j). As an example, 
here we expand 1(1 + e) f|4.2j) up to e 5 ; the result agrees with (|4.8p . (|6.7p . The 
following text is a Mathematica notebook exported to ETgX and slightly edited for 
readability. 



Initialization 

« HypExp 

*_*_*_*_*_* HPL 2 *_*_*_*_*_* 

Author: Daniel Maitre, University of Zurich 

Rules for minimal set loaded for weights: 2, 3, 4, 5, 6, 7, 8. 

Rules for minimal set for + - weights loaded for weights: 2, 3, 4, 5, 6, 7, 8. 

Table of MZVs loaded up to weight 8 

Table of values at I loaded up to weight 7 

SHPLFunctions gives a list of the functions of the package. 

SHPLOptions gives a list of the options of the package. 

More info in |hep-ph/0507152[ |hep-ph/0703052| and at 

http:/ /krone. physik.unizh.ch/~maitreda/HPL/ 

*********** HypExp 2 ************ 
Authors: 

Tobias Huber: RWTH Aachen, 

Daniel Maitre: SLAC, University of Zurich. 

HypExp loaded! It allows the expansion of hypergeometric functions around their parameters. 
The new provided commands are: 
- HypExp 
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- HypExpInt 

- HypExpU 

- HypExpAddToLib 

- HypExpIsKnownToOrder 

More info in |hep-ph/0507094| and at 

http: / /krone. physik.unizh.ch/~maitreda/HypExp/ 

F = Gamma; 



Kotikov 

We start from various expressions for I(n + e) (15. 9p with n = 0, 1, 2 as linear 
expressions in a hypergeometric function F, and divide them by 

r(l + 3e)r 2 (l-e)r(l-2g) _ 
r(l + e)r(l-4e) ' 

then their expansions don't contain the Euler constant 7, and the right-hand sides of 
the IBP relations (|2.2p become rational functions. Then we use these IBP relations 
to reduce the cases n = 0, 2 to n = 1, and multiply by the appropriate factor to 
obtain I (|4,2p . The result contains 

$80. = C53 + Cs • 

x_, f_, n_]:=Module[{a, b, y}, a = D[x, F];b = x/.F — > 0; b = b/a; 
a = a * T[l + e] * P[l - 4 * e]/(T[l + 3 * e] * T[l - e] A 2 * T[l - 2 * e]); 
y = Series [HypExpf/, e, n], {e, 0, n}]; y = MapfExpand, FunctionExpand[a *(y + b)]]; 
y = Switch[j, 0, 

Map[Expand, ((1 - 3 * e) * y - (1 - 5 * e)/(3 * e A 2 * (1 - 3 * e) * (1 - 4 * e)))/(2 * e)], 
2, Map [Expand, (-(1 + 2 * e) * y + (1 + 5 * e)/(3 * e A 2 * (1 + e) * (1 + 2 * e)))/(3 * e)], 

-»»]; 

Map[Expand, FunctionExpand[y * (1 — 2 * e) * T[l — e] A 2* 
Sqrt[r[l - 2 * e] * T[l + 3 * e]/(T[l + e] * T[l - 4 * e])]]]] 

FK[n_] = HypergeometricPFQ [{ 1 , 2 - 2 * e, n + 2 * e}, {n + 1 + e, n + 1 + 2 * e}, 1]; 

j£T[n_] = 2 * r[l - e] * T[l - n - 2 * e] * T[n - 1 + 3 * e]* 

(it * Cot[7r * (n + 3 * e)]/r[2 - 2 * e]- 

P[l - e]/((n + 2 * e) * T[2 - n - 4 * e) * T[n + 1 + e]) * F); 

Il = i[l,tf[l],FK[l],7] 

6C(3) + ^ + 102 C (5) £ 2 + - 24 C (3) 2 ) + (l413 C (7) - ^® ) + 

/648HPLs8a , v.,* 6017vr 8 \ - _ , fi . 

(— 228<(3)C(5) + — ) £ 5 + 0( £ 6 ) 

{x[0,K[0],FK[0],8] -Il,i[2,K[2],FK[2],6] -11} 

{0(^,0^)} 
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Broadhurst, Gracey, Kreimer 

Here we use (|6.8p - (l6.1ip . In two cases, I was unable to expand up to e 5 : after running 
for a long time, Mathematica said "No more memory available". In these cases, I 
expanded up to e 4 . 

6[n_]:= - 2 * T[l - e] A 2 * T[n - 1 + 3 * e]* 

r[l - n - 2 * e]/(e * (1 - 2 * e) * T[n + e] * T[2 - n - 4 * e\) 

a[n_]:=r[n + e] * T[2 - n - 4 * e]/(r[l - e] * T[l - 2 * e]) 

Fl[n_]:=HypergeometricPFQ[{l, — e, n — 1 + 3 * e}, {n + e, n + 3 * e}, 1] 

Bl[n_]:=6[n] * ((1 — n — 4 * e)/(l — n — 3*e) * F — a[n] *ir*e* Cot[3 * ir * e) — 1) 

{i[l,Bl[l]-Il,Fl[l],8],i[0,Bl[0],Fl[0],6]-Il,i[2,Bl[2],Fl[2],7]-Il} 

{0(£ 6 ),0(, 5 ),0( £ 6 )} 

F2[n_]:=HypergeometricPFQ[{l, 1 - n - e, 1 - n - 3 * e}, {1 + e, 2 - n - 3 * e}, 1] 

B2[n_]:=6[n] * (—(1 — n — 4 * e)/(l — n — 3*e)*F + a[n] * 7r * e * Cot[7r * e] — e/(l — n — 3 * e)) 

{i[l,B2[l],F2[l],8] -Il,i[0,B2[0],F2[0],7] - Il,i[2,B2[2],F2[2],7] -11} 

{0(e 6 ),0(, 6 ),0(, 6 )} 

F3 [n_] :=HypergeometricPFQ [{ 1 , n - 1 + 4 * e, n - 1 + 2 * e}, {1 + e, n + 2 * e}, 1] 
B3[n_]:=6[n] * (— (n — 1 + e)/(n — 1 + 2 * e) * F + a[n] *ir*e* Cot[7r *e]+ e/(l — n — 2 * e)) 

{*[1,B3[1],F3[1],8] -I1,»[0,B3[0],F3[0],7] - Il,i[2,B3[2],F3[2],7] -11} 

{0( £ 6 ),0(, 6 ),0( £ 6 )} 

F4[n_]:=HypergeometricPFQ[{l, -e, 1 - n - 2 * e}, {2 - n - 4 * e, 2 - n - 2 * e}, 1] 
B4[n_]:=6[n] * ((n — 1 + e)/(n — l + 2*e)*F + a[n] * ir * e * Cot[2 * ir * e] — 1) 
{i[l,B4[l],F4[l],8] - Il,i[0,B4[0],F4[0],7] -I1,»[2,B4[2],F4[2],6] -11} 

{0( £ 6 ),0(, 6 ),0(, 5 )} 

Hathrell 

For this particular problem we can also use (IA.2|) . 

FH = HypergeometriePFQ [{ 1 , 2 - 2 * e, 2 - 3 * e}, {3 - 4 * e, 3 - 3 * e}, 1]; 

H = -2 * e * T[l - e] * T[-l + 3 * e] * T[l - 2 * e]/T[l + e]* 

(T[l - e]/((2 - 3 * e) * T[3 - 4 * e]) * F - T[-l + 2 * e] * T[e] * Cos[2 * tt * e]); 

t[0,fT,FH,6] -II 
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Kazakov 



Expanding hyper geometric functions of —1 is a little tricky. First, we should switch 
off automatic conversion of harmonic polylogarithms to usual polylogarithms (oth- 
erwise values of polylogarithms on their cuts will be generated). These expansions 
involve, in addition to multiple £ values, alternating (Euler-Zagier) sums. After 
reducing them to a minimal set, many new constants appear (including log 2 and 
polylogarithms of |). These new constants cancel in the square bracket in (|3.6p (it 
is calculated by the function br). 

$HPLAutoConvertToKnownFunctions = False; 

FK 1 [n_] :=HypergeometricPFQ [{ 1 , 2 - 2 * e, n + 2 * e}, {1 + e, n + 1 + 2 * e}, -1] 
FK2[n_]:=HypergeometricPFQ[{l, 2 - 2 * e, 2 - n - 3 * e}, {1 + e, 3 - n - 3 * e}, -1] 



brp_, n_]:=Module[{fl, £2, nl, n2, y}, nl = n - If[j===2, 1, 0]; 
n2 = n-If(j ===0,1,0]; 

fl = Map [Expand, Series[HypExp[FKl[?], e, nl], {e, 0, nl}]]; 
f2 = Map[Expand, Series [HypExp[FK2 [?'], e, n2], {e, 0, n2}]]; 
Map [Expand, fl/(j + 2 * e) + f2/(2 - j - 3 * e)]] 

ik[j_, n_]:=Module[{y}, y = br\j, n]; 
y = Map [Expand, 

FunctionExpand[y - T[l + e] * T\j + s] * T[2 - j - 4 * e]/T[2 - 2 * e] * Cos[tt * e]]]; 
y = Map [Expand, 2 * Pochhammer[l + 3 * e, j — 2]* 

Pochhammerfl — 2 * e, — j]/ (Pochhammer[e, j] * Pochhammerfl — 4 * e, 1 — j]) * y]; 
y = Switch [7, 0, 

Map[Expand, ((1 - 3 * e) * y - (1 - 5 * e)/(3 * e A 2 * (1 - 3 * e) * (1 - 4 * e)))/(2 * e)], 
2, Map [Expand, (-(1 + 2 * e) * y + (1 + 5 * e)/(3 * e A 2 * (1 + e) * (1 + 2 * e)))/(3 * e)j, 
-tV]; 

Map[Expand, FunctionExpand [y * (1 — 2 * e) * T[l — e] A 2* 
Sqrt[r[l - 2 * e] * T[l + 3 * e]/(T[l + e\ * F[l - 4 * e])]]]] 

{ik[l, 8] - II, ik[0, 8] - II, ik[2, 8] - 11} 

(0(e 6 ),0(e 6 ),0(e 6 )} 
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